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Abstract 

We study arrival directions of 1.4 ■ 10^ extensive air showers (EAS) registered with the EAS-1000 
04 I Prototype Array in the energy range 0.1-10 PeV. By applying an iterative algorithm that provides 

uniform distribution of the data with respect to sidereal time and azimuthal angles, we find a 
I number of zones with excessive flux of cosmic rays (CRs) at > 3a level. We compare locations of 

^ ■ the zones with positions of galactic supernova remnants (SNRs), pulsars, open star clusters, OB- 

OO I associations, and regions of ionized hydrogen and find remarkable coincidences, which may witness 

CO i in favour of the hypothesis that certain objects of these types, including the SNRs Cassiopeia A, 

the Crab Nebula, the Monogem Ring and some other, provide a noticeable contribution to the flux 
of CRs in the PeV range of energies. In addition, we flnd certain signs of a contribution from the 
I M33 galaxy and a number of comparatively nearby groups of active galactic nuclei and interacting 

■ galaxies, in particular those in the Virgo cluster of galaxies. The results also provide some hints 

I for a search of possible sources of ultra-high energy (UHE) cosmic rays and support an earlier idea 

that a part of both UHE and PeV CRs may originate from the same astrophysical objects. 

o 

to '. TMTOWTDI* 
^ ' The Perl community mantra 



1. Introduction 

In spite of the fact that cosmic rays were discovered more than 90 years ago [1] , the problem of their origin for 
' energies greater than 100 TeV remains unsolved. One of the important directions in numerous approaches 
to the problem is an analysis of arrival directions of CRs. Such an analysis has been performed with the 
data sets obtained with practically all extensive air shower experiments, see [2], [3], [4], and [5] for the latest 
reports and a list of references. One of the main results of the majority of these investigations is that there 
is no signiflcant anisotropy in the energy range 0.1-100 PeV [6]. 

Another intriguing and long-standing problem in cosmic ray physics is the so called 'knee' around 
3 • 10^^ eV in the CR energy spectrum [7]. The knee is a point where the spectral index of the all-particle 
differential power-law spectrum changes from approximately —2.7 to —3.1. There are a number of models 
aimed to explain this feature (as well as the origin of CRs), see, e.g., [8], [9], [10], and [11] for a number of 
recent works and [12] for an earlier review. Still, none of them seems to be fully established yet. One of the 
modern 'astrophysical' approaches to the problem is the 'single source model' by Erlykin and Wolfendale [13]. 
The model explains the knee as a result of a contribution from one comparatively recent (~ 100 kyr) and 
nearby (~ 300 pc) supernova remnant accompanied by a pulsar. Evidently, an analysis of arrival directions 
of CRs in the energy range around the knee can be useful as a (partial) experimental test of the model. 

An interest to the single source model has recently received an additional impulse after the discovery 
that the pulsar PSR B0656-I-14 is located near the center of the Monogem Ring SNR at the distance of 
about 290 pc [14]. Since then, there has been a number of reports on the analysis of the flux of cosmic 
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rays with energies around the knee from the direction to this source. Interestingly, some of the results are 
controversial. Namely, Chilingarian et al. [2] and Benko et al. [4] find a significant excess of EAS from zones 
near the pulsar PSR B0656+14 while the KASCADE collaboration considers this excess as negligible [3]. To 
shed an additional light on the situation, we perform an analysis of arrival directions of EAS registered with 
the EAS-1000 Prototype Array, which has been operating at Moscow State University. Unfortunately, the 
energy range covered by the array is not sufficient for a comprehensive investigation of arrival directions of 
EAS around the knee since it mostly covers a range of energies just below the knee. Still, we find the first 
results of this work sufficiently interesting and promising to report. 

The most important result is the discovery of a whole number of zones with an excessive flux (ZEE) of 
CRs around 1 PeV at > 3a level. The majority of the ZEE correspond to the locations of galactic SNRs, 
pulsars, open star clusters, OB-associations, and regions of ionized hydrogen. In particular, we confirm the 
conclusions of Chilingarian et al. [2] and Benko et al. [4] on the excessive flux of CRs from the direction to 
the Monogcm Ring. Still, a number of the zones have no or just a few objects of the above types nearby. 
We find that one of these 'empty' or 'underfilled' ZEE contains the M33 galaxy inside, while the majority 
of others have neighbouring active galactic nuclei (AGN) and/or interacting galaxies at redshifts z < 0.01 
with a big group of them located in the Virgo cluster of galaxies. In addition, we find a few remarkable 
coincidences between positions of the ZEE and arrival directions of ultra-high energy cosmic rays registered 
with the AGASA array [15]. 

2. Experimental Data 

The EAS 1000 Prototype Array consists of eight detectors situated in the central part of the EAS MSU 
array along longer sides of the 64 m x 22 m rectangle [16]. The array is located at 37°32.5'E, 55°41.9'N at 
approximately 200 m above sea level. 

The data set under consideration includes 1,668,489 EAS registered during 203 days of operation of the 
array in the period from August 30, 1997, till February 1, 1999. The arrival directions are determined for 
1,366,010 EAS. A number of charged particles (electrons) is found for 826,921 EAS. It happens that A^e 
of 95.2% of showers with zenith angles 6 < 45° lie in the range 10^-1.1 • 10^ particles with iVc = 1.2 • 10^ and 
the median value equal to 6.0 • 10^. Thus we estimate that the overwhelming majority of primary cosmic 
rays that give life to the EAS in the data set have energies in the range E » 0.1-10 PeV. Still, only 92,212 
EAS have A^e > 3 • 10^. Thus one may treat the data set as mostly covering an interval of energies just below 
the knee. 

3. Method of Data Analysis 

Perhaps the majority of investigations of anisotropy rely on the Rayleigh method for the calculation of the 
amplitude and the phase of the first and (sometimes) the second harmonics [6] . The method by itself is based 
on the Fourier transform. This does not seem to be a fruitful approach in all situations. For example, one 
can easily check that an analysis of the one-dimensional Fourier spectrum does not necessarily give reliable 
results if the amplitude of a periodic signal is less than or of the order of 0.3% of that of noise for a sample 
that consists of 10^ points. Thus we expect that a similar situation can take place in the case of anisotropy 
analysis. This makes us develop another approach, which may be considered as an advanced version of the 
method employed in [17]. The algorithm consists of three main steps: data alignment, data averaging, and 
selection of ZEE. 

Step 1: Data Alignment. The main aim of the step is to obtain a data set that consists of EAS with 
a uniform distribution of azimuthal angles (f) and sidereal time s. 

1. Fix values Acf) and A.s of the width of bins in histograms of 4> and s. Produce a histogram of </>. By Ni 
denote the number of EAS in the ith bin. 

2. In the histogram of 0, find a bin with the least number of showers, Nj^in- At random, exclude (A^i — A^min) 
showers from each bin. 

3. Obtain a histogram of s for the remaining data set. 'With this histogram, perform an 'alignment' of 
data exactly as in item 2. 

4. Produce a histogram of ^ for the data set obtained in item 3. Finish the procedure if 

max lA^j - iVl < 2.5, and max |A^j - iVl < 3a, (1) 

i i 
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where N is the mean number of EAS in the bins of the histogram, and a is the standard deviation. 
Otherwise, repeat items 2 and 3. 

Conditions (1) were chosen after comparing a number of empirical criteria. In practice, the second 
inequaHty is always fulfilled earlier than the first one and may be omitted. 

We stress that the data set obtained as a rcsiilt of Step 1 contains equal number of EAS in each bin of 
the histogram of the sidereal time s thus providing uniform distribution of the data with respect to s at the 
chosen scale As. The same is approximately true for the distribution with respect to (f>. 

For the data set under consideration. Step 1 leads to an exclusion of up to 20% of all EAS in the data 
set. It is obvious that results of an analysis of the 'aligned' data set will strongly depend on the choice of 
the excluded EAS. Hence we come to the necessity to implement Step 1 numerous times thus obtaining a 
family of 'aligned' data sets. Each time, an independent sample of EAS is excluded from the whole data set. 

Step 2: Data Averaging. For each 'aligned' data set, produce a map of the distribution of EAS 

arrival directions in equatorial coordinates. Each map consists of 1° x 1° cells ('unit' cells) and covers the 
complete observable region with the right ascension a = ... 24 h and the declination S > —13°. For each 
cell in all maps, a number of EAS inside the cell is found. After this, an averaged map is produced. 

At this point, we arc ready to join unit cells into larger ones and look for cells with an excessive flux 
(CEF) of cosmic rays. We do not analyse unit cells because it is seems unlikely that charged particles of 
PeV energies preserve directions to their sources with such precision. 

Step 3: Selection of CEF. 

1. Divide the averaged map into strips of equal width AS > 3° in the declination. The boundaries of strips 
of equal width are shifted with respect to each other by 1°. 

2. Split each strip into colls of equal width Aa. For each strip and each partition into cells, find the mean 
number N of EAS for all cells within the strip and the standard deviation a. Mark a cell with N* EAS 
inside as having an excessive CR flux if N* — N > 3.1a, and A''* > 10. 

Remark 1. It often happens that the algorithm selects a number of intersecting or overlapping CEF, so 
that it is inconvenient to discuss each of them separately. Instead, we shall discuss zones of excessive flux 
(ZEE), by which wc mean sets of intersecting or overlapping CEF. For simplicity, we shall use this term even 
in case a ZEE consists of a single cell. 

Remark 2. It is not important for the algorithm how the values of Aa and A^ are chosen. The majority 
of ZEF discussed below are combined of cells such that Aa equals a rounded value of A(5 / cos (5, where 5 is 
the mean value of the declination for the strip. This guarantees that cells with the same A(5 but located at 
different declination, have an approximately equal area. This is a traditional choice for similar investigations. 
In what follows, we call these cells 'regular'. Still, there are no reasons to restrict the analysis to regular 
cells only. Thus, we have also studied cells with arbitrary values of Aa providing that they are integer 
dividers of 360°. These 'irregular' cells are in no sense worse than 'regular' ones. More than this, in some 
cases an employment of irregular cells allows one to find CEF that are more pronounced than their regular 
counterparts in the sense of deviation from N . On the other hand, for any regular CEF, there is an irregular 
one with the same or close size. A number of interesting irregular CEF will be presented below. 

It is important to mention that all CEF have a neighbourhood with a flux of CRs that is higher than 
the average for the strip though it does not reach the level of S.lcr. In our opinion, this means that CEF are 
not just random fluctuations of the level of the flux caused by the algorithm. 

Remark 3. The rule used to choose CEF implicitly assumes that the number of EAS within cells of 

any strip obey the Gaussian distribution with the corresponding parameters. The hypothesis seems to be 
possible since these numbers are generically non-integer due to the averaging of the data set. Nevertheless, 
one may choose another rule to select CEF. 

Remark 4- An analysis of the algorithm described above has revealed that maps of 'averaged' EAS 
arrival directions obtained with 100 cycles of alignment differ very little from each other. Still, in order to 
make our conclusions even more robust, we present results obtained for 1000 cycles. To perform calculations, 
we have employed GNU Octave [18] running in Linux. 

Remark 5. In what follows, we present results obtained for A(j) = 20° and As = 3 h, which provide 
a comparatively small number of excluded EAS and an acceptable speed of the algorithm. An analysis of 
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averaged maps obtained for different values of Acj) and As in the range = 10° . . . 30° and As = 1 ... 3 h 
and cells of the size 5° x 5° has revealed that the results differ from those obtained for A(p = 20° and 
As = 3 h by at most O.lcr. Thus one may expect that the condition N* — N > 3.1a used to select CEF 
guarantees that N* — N > 3.0a for other sensible values of Acf), As, and sizes of cells. 
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Figure 1: Zones with excessive flux of cosmic rays (rectangles) and galactic SNRs (o) in the region 
S > —15°. The ZEF are numbered from top to bottom according to the values of their upper 
boundaries. The curves show the galactic plane. Filled circles show SNRs located within 2°-neigh- 
bourhoods of the ZEF. 



4. Analysis of the Data Set 

The above procedure applied to the whole data set resulted in the 'averaged' set of 1,134,325 EAS and 193 
regular and an even greater number of irregular CEF. The regular CEF and a number of irregular CEF were 
joined into 37 ZEF, see Fig. 1. A number of CEF are embedded into each other or shifted with respect to 
each other by 1°. To make this and the following figures as clear as possible, we omit some of these CEF 
but show their joint boundaries instead. 

An important question is how many EAS are there inside ZEF? Below we shall give exact numbers but 
it may be useful to remember the following: (i) any CEF contains at least 10 EAS; (ii) any unit cell located 
in the region 26° < S < 80° contains > 10 EAS; (iii) the mean number of EAS inside unit cells located in 
the region 19° < 6 < 85° is greater than 9; (iv) min^ = -12.38° but there are only 50 EAS with S < 0°. 

4.1. Galactic SNRs. Let us begin an analysis of positions of possible CR sources with galactic supernova 
remnants. For this purpose, we use the January 2004 version of the D. Green Catalogue [19]. 

Figure 1 shows positions of 96 galactic SNRs located in the region S > —15°. (The Monogem Ring 
remnant is not shown since it is not present in the main list of [19].) It happens that 15 SNRs lie inside or 
at angular distances A < 2° from the nearest ZEF. Namely, 

• ZEF No. 5 has the SNR G132.7+1.3 (HB3) inside and the remnant G130. 7+3.1 (SN1181) located at 
the angular distance A — 1.1°; 

• ZEF No. 6 has the SNR G127.1+0.5 (R5) at A = 0.9°, the SNR G126.2+1.6 at A = 2.0°, and the 
remnants associated with ZEF No. 5 at A = 1.0° and A = 1.8° respectively; 

• ZEF No. 8 has the SNRs G109.1-1.0 (CTB 109) inside, G106.3+2.7 at A = 0.8°, and Gill. 7-2.1 
(Cassiopeia A) at the angular distance A = 2.0°; 

• ZEF No. 18 has five SNRs inside: G179.0+2.6, G180.0-1.7 (S147), G182.4+4.3, G184.6-5.8 (the Crab 
Nebula, SN1054), and G189. 1+3.0 (IC443); 
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• ZEF No. 20 has the SNR G65.3+5.7, located at A = 1.2° and the SNR G59.8+1.2 at A = 1.8°. 

• ZEF No. 35 has the SNR G206.9+2.3 (PKS 0646+06), located at A = 0.8°. 



Table 1. Some parameters of the ZEF with galactic SNRs inside or at angular distances A < 2°. 
Notation by columns: 'ZEF' is the munbcr of the ZEF as shown in Fig. 1 (letters 'Ir' mark irregular 
CEF), a° and 5° are the ranges of values of the right ascension and the declination for the listed cells 
respectively, N* is the number of EAS in the cell, N is the mean number of EAS in the corresponding 
strip for the given value of Aa, a is the standard deviation for the strip, D = {N* — N)/a, i.e., the 
deviation of N* from N in units of a. 
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Notice that a number of these SNRs are widely discussed in literature as possible sources of CRs in the 
region of the knee. Among them, Cassiopeia A (see, e.g., [20], [21] and references therein), the Crab Nebula, 
and the SNRs S147 and G65.3+5.7 [22], [23]. 

An extension of the selection zone up to the angular distance A < 3° adds two more SNRs to the list, 
namely G59. 5+0.1 and G63. 7+1.1 both lying near ZEF No. 20. Another prominent SNR, the Monoceros 
Nebula, is located at A = 3.2° to the west from ZEF No. 35. 

Table 1 contains some parameters of the discussed ZEF. Among them, ZEF No. 5, 6, and 20 are made 
of irregular cells. In the boundaries shown in Fig. 1, ZEF No. 6, 8 and 35 consist of 18, 10, and 4 CEF 
respectively. For the sake of brevity, only the most interesting of them are presented in Table 1. Notice 
that all these ZEF except No. 35 lie in the region with a high density of EAS per unit cell and thus contain 
hundreds and even thousands of showers. 

4.2. Pulsars. Figure 2 shows the mutual position of the ZEF and pulsars known in the region S > —15° [24], 
version of October 11, 2004. Forty-six of these 640 pulsars lie inside the ZEF, other 34 pulsars lie at angular 
distances A < 2° from the boimdaries of the nearest ZEF. As one can see from the figure, the list of zones 
with 'neighbouring' pulsars consists of 21 ZEF: No. 2-6, 8, 9, 14, 16, 18, 20, 21, 26, and 29-36. Notice that 
all ZEF discussed above in connection with galactic SNRs belong to the list. Some parameters of the ZEF 
'selected' by pulsars are presented in Tables 2 and 3. The biggest number of pulsars within 2°-neighbourhoods 
are found for ZEF No. 32 (26 PSRs), No. 8 and 18 (8 PSRs), No. 6 and No. 20 (7 and 6 PSRs respectively). 
Each of ZEF No. 5, 9, 14, and 21 has four neighbouring pulsars. It is interesting to mention that 13 of 640 
pulsars in the region 6 > 15° are associated with galactic SNRs. Five of them lie within A < 2° from the 
ZEF. 

We have already mentioned the pulsar PSR B0656+14 in connection with the single-source model by 
Erlykin and Wolfendale [13]. The pulsar has equatorial coordinates a = 104.95°, S = 14.24° and is located 
close in projection to the center of the Monogem Ring SNR [14]. As is clear from Fig. 2 and Table 3, the 
pulsar lies at the angular distance A = 2.05° from the left boundary of the outer cell of ZEF No. 29. This 
is an irregular cell of the size 9° x 4°. The interior (regular) 3° x 3° cell is located at A = 7.05° from the 
pulsar. 

The interior cell is especially interesting in connection with the article by Chilingarian et al. [2], in 

which arrival directions of more than 2 ■ 10^ EAS with the size Nc > 10^ detected by the MAKET-ANl 
experiment have been analysed. As a result, a 3° x 3° CEF located at a = 111° . . . 114°, 6 = 12.5° . . . 15.5° 
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Figure 2: The ZEF and pulsars in the region 6 > —15° [24]. Different symbols denote pulsars 
located at different distances d from the solar system: o — d < 3 kpc, A — 3 < < 6 kpc, x — 
d > 6 kpc or the distance is unknown. Filled circles, filled triangles, and asterisks (*) are used 
respectively to show pulsars located within 2°-neighbour hoods of the ZEF. 

has been found. Evidently, this CEF is very close to the interior cell in ZEF No. 29. Here one should take 
into account that (i) we have not studied cells with boundaries located at non-integer values of a or S, and 
(ii) Chilingarian ct al. have not studied cells with the size larger than 3° x 3°. Thus we think this is a 
remarkable coincidence of the results, obtained in different experiments and by different methods of data 
analysis. 

It is worth mentioning that the Monogem Ring remnant covers a huge region with the diameter of 
approximately 25° with its most 'pronounced' part being at the right hand side of the circle, see the figure 
in [14]. (To the contrary to our figures, Thorsett et al. [14] use the grid with bigger values of a located at the 
left-hand side of the figure.) Thus it is clear that ZEF No. 24 and 35 are intersected by the most pronounced 
part of the Monogem Ring, cf. Fig. 1. 

As one can see in Fig. 2, two other pulsars are located comparatively close to ZEF No. 29. These are 
J071H-0931, located on the top of ZEF No. 35, {d = 2.39 kpc, A = 2.5° from ZEF No. 29) and J0751+1807 
{d = 2.02 kpc, A = 2.7°), which fies to the NE from ZEF No. 29. (This pulsar looks to be located closer to 
ZEF No. 24 but in fact the corresponding angular distance equals nearly 2.8°.) 

Another interesting issue is the distribution of distances d to the 'neighbouring' pulsars. As we have 
already mentioned above, there are 80 pulsars that belong to the 2°-neighbourhoods of the ZEF. One of 
them lies at c? > 10 kpc. Distances to four other pulsars are unknown. The remaining 75 pulsars are located 
at d < 10 kpc. As for the whole set of 640 pulsars in the region 5 > —15°, 582 of them have c? < 10 kpc, 
and d > 10 kpc for 49 pulsars. Distances to the remaining 9 pulsars are unknown. Thus let us take a look 
at the distribution of distances < 10 kpc. To do this, define 



K{d) 



582 N^^°^{d) 



(2) 



where Npg^ and -/Vpgj^ are the number of pulsars within the corresponding intervals of d for the 2°- 
neighbourhoods of ZEF and for all pulsars respectively. A value K = 1 corresponds to the case when a 
fraction of neighbouring pulsars at the given interval of d equals that for the whole set of pulsars. 

Figure 3 shows the behavior of K{d) in 0.5-kpc intervals. It is clearly seen that a fraction of neighbouring 
pulsars with d > 8 kpc is negligibly small. On the other hand, K > 1 for d < 4 kpc with the most pronounced 



6 



Table 2. Some parameters of the ZEF with pulsars inside. For composite ZEF made of multiple 

CEF, the corresponding number is given in parenthesis after the number of the zone. See Table 1 
for other notation and parameters of ZEF No. 5, 6, 8, 18, and 20. 
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bins located at 1.5 < d < 2 kpc and d < 0.5 kpc.^ In our opinion, the latter fact can witness in favour of the 
hypothesis that a noticeable fraction of the overall flux of CRs near the knee is produced by nearby pulsars. 

Totally, we find five pulsars located at distances d < 0.5 kpc from the solar system that belong to the 
2°-neighbourhoods of ZEF: 

• PSR B0656+14 at A « 2° from ZEF No. 29 {d = 0.29 kpc); 

• PSR .11136+1551 inside ZEF No. 30 {d = 0.36 kpc), also at A = 1.5° from ZEF No. 14; 

• PSR J1744-1134 at A = 1.6° from ZEF No. 32 {d = 0.36 kpc); 

• PSR J0814+7429 inside ZEF No. 2{d = 0.43 kpc); 



^ The ATNF database [24] gives two 'basic' values of distance to a pulsar: the best estimate of the pulsar 
distance and the distance based on the Taylor & Cordes electron density model [25]. For definiteness, we 
use the best estimate. Nevertheless, our results do not change qualitatively if the second value is used. In 
this case, the highest bins correspond to distances d < 0.5 kpc and 1.5 < d < 2.5 kpc. 
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Table 3. Some parameters of the ZEF with pulsars at angular distances < A < 3°. See Table 1 
for notation and parameters of ZEF No. 35. 
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Figure 3: The behaviour of the function K{d), see Eq. (2), for 0.5-kpc intervals. 



• PSR J0700+6418 inside ZEF No. 3 {d = 0.48 kpc). 

The Geminga pulsar, which is often discussed as a possible source of PeV CRs {d = 0.16 kpc) is located 
at the angular distance A = 3.9° from the bottom-right corner of ZEF No. 18. 

Rem,ark 6. Wc do not really need irregular CEF for zones No. 4 and 33 to obtain 80 pulsars in the 
2°-neighbourhoods of the ZEF. If we exclude them from consideration, the corresponding pulsars 'move' 
from inside the ZEF to A = 1.7° and A = 0.8° for ZEF No. 4 and 33 respectively. 

As we have stated above, there arc 21 ZEF that have pulsars within 2°-neighbourhoods. Three more 
zones (No. 17, 24, and 27, cf. Fig. 1) are added to the list if we select pulsars located at A < 3°. (In this 
case, we find 23 more neighbouring pulsars.) An extension of the selection region up to A < 4° does not 
change the list of 'selected' ZEF. Hence, there remain 13 ZEF without neighbouring pulsars and/or galactic 
SNRs. It is not a simple task to estimate a possible excess of the flux of CRs due to a contribution from a 
pulsar, see, e.g., [26], but the number of neighbouring pulsars for some of the 'selected' ZEF does not look to 
be 'sufficient' either. In particular, this applies to huge zones No. 14, 16, and 30 with the first two of them 
lying in the region with a comparatively big number of EAS per unit cell, see Table 2. Thus we are lead to 
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the necessity to analyse positions of other possible sources of PeV CRs. 




4.3. Open Star Clusters and OB- associations. The fact that most open star clusters (OSCs) are dominated 
by the hottest, i.e., 0-type and B-type, giant blue stars made us take a look at positions of these objects. 
We have analysed coordinates of 299 OSCs that lie in the region 6 > —15° and seem to be located in the 
Milky Way Galaxy. The list of these OSCs was mainly prepared with the use of the SIMBAD database [27] . 
The result is shown in Fig. 4: 42 OSCs lie inside the ZEF, 32 more clusters are located at angular distances 
< A < 2° from the ZEF. 

Similar to the cases of pulsars and galactic SNRs discussed above, ZEF No. 5, 6, 8, and 18 are among 
the 'leaders' with the biggest number of objects lying within 2°-neighbourhoods. What seems to be more 
interesting is that the hst of 'selected' ZEF is now enriched with ZEF No. 12 (OSCs C 0417+448 and 
NGC 1605), ZEF No. 13 (OSC Czernik 35), and ZEF No. 25 (NGC 1647). In addition, ZEF No. 24, which 
has the J0751+1807 pulsar at the angular distance A = 2.8°, now gets an object inside (NGC 2420). ZEF 
No. 15 is added to the list if we extend the selection region up to A = 3°, cf. Fig. 1. In this case, ZEF No. 25 
and 26 also obtain a neighbouring OSC, namely the Hyades (C 0424+157). It is worth mentioning that the 
Hyades is a nearby OSC, located at the distance of about d = 46 pc from the solar system. 

It is interesting that ZEF No. 24, 29, and 35, discussed above in connection with the Monogem Ring 
SNR, contain 1, 2, and 3 OSCs respectively inside. A small cell at the NW corner of ZEF No. 18 also 
becomes 'selected' for the first time. In particular, it contains NGC 1931, M36, and M38 inside and M37 
just nearby. (These three Messier objects are located at d = 1.26-1.35 kpc.) Finally, let us mention an OSC 
that lies near the center of ZEF No. 14. This is the Coma star cluster, located at d = 90 ± 5 pc from the 
solar system. 

Table 4. Parameters of ZEF 'selected' by OSCs. See Table 1 for notation. 
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Thus an assumption that open star clusters can be sources of PeV CRs adds three or four more ZEF to 
the hst of zones with neighbouring CR sources. Their parameters are hsted in Table 4. 

The so-called superbubble model [28] stimulated us to analyse positions of OB-associations. We have 
found that at least seven OB-associations are located at angular distances A < 2° from ZEF No. 6, 8, and 36. 
Three other associations are located at A « 2.2°-2.5° from ZEF No. 20. Remarkably, open star clusters and 
OB-associations taken together with galactic SNRs 'select' all ZEF in the vicinity of the galactic plane. 




20 40 60 80 100 120 140 160 180 200 220 240 260 280 300 320 340 360 



Figure 5: The ZEF and galactic HII regions (+) [29]. Filled circles show regions lying within 
2°-neighbourhoods of the ZEF. 



4.4. Galactic HII Regions. Another group of astrophysical objects that may be considered as possible 
sources/ accelerators of cosmic rays consists of regions of ionized hydrogen (HII). We have analysed positions 
of galactic HII regions given in the catalogue by Paladini et al. [29]. The result is shown in Fig. 5. 

Seventy-three of 728 HII regions with declination S > —15° belong to the 2°-neighbourhoods of the 
ZEF. Thirty-seven of these regions lie inside the corresponding ZEF. The biggest number of HII regions is 
found in the vicinity of ZEF No. 8 (29 regions, 17 of them inside the zone), ZEF No. 18 (15 and 11 regions 
respectively), and ZEF No. 5 (9 and 2 regions respectively). Similar to the case of open star clusters, ZEF 
No. 12 and 25, which do not have any galactic SNRs or pulsars in their vicinity, become 'selected' with 
respectively three and one HII regions inside. An additional 'support' is also provided to ZEF No. 34 (3 HII 
regions nearby) and to the small cell in the NW corner of ZEF No. 18. Still, the analysis of positions of 
galactic HII regions has not added new zones to the list of ZEF with possible sources and/or accelerators of 
PeV CRs in their vicinity. 

An analysis of positions of Wolf-Rayet (WR) stars did not reveal anything remarkable except that there 
are more than 100 WR stars within the boundaries of a small 5° x 4° cell located inside ZEF No. 16. We 
shall discuss this fact later. At the moment, we see that 25 of 37 ZEF have at least one galactic object of the 
types discussed above in their 2°-neighbourhoods, 28 ZEF if we consider 3°-neighbourhoods. Thus, there 
remain at least 9 ZEF (No. 1, 7, 10, 11, 19, 22, 23, 28, and 37) that do not have known galactic objects of the 
above types in their vicinity. A number of the zones may look 'underfilled'. If we believe that presented ZEF 
somehow reflect a situation with the angular distribution of the flux of PeV cosmic rays then the appearance 
of 'empty' and 'underfilled' ZEF raises a question on other possible CR sources. An unexpected hint comes 
from an analysis of arrival directions of the ultra-high energy cosmic rays (UHECRs) registered with the 
AGASA array. 

Earlier, we have used the same data set for a study of arrival times of EAS. In particular, we found 20 
groups ('cluster events') of consecutive EAS that produced bursts of the count rate [30], [31]. The majority 
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a° 

Figure 6: The ZEF and arrival directions of the UHECRs with E ^ 4 ■ 10^^ eV registered with the 
AGASA array (+) [15]. FiUed circles mark UHECRs lying within 4°-neighbourhoods of the ZEF. 



of EAS within these groups did not have close arrival directions but certain clusters covered approximately 
the same regions of the sky. Surprisingly, we have found a group of four such clusters with the famous 
AGASA 'C2' triplet near the center of the corresponding region, see Fig. 30 in [31]. This gave rise to an 
idea that there may be some connection between UHECRs and cosmic rays in the PeV range. The present 
investigation gave us a chance to test the idea. 

Figure 6 shows arrival direction of 58 UHECRs with energies E > 4-10^^ eV registered with the AGASA 
array [15]. Twenty-one of these UHECRs belong to 4°-neighbourhoods of the ZEF. It is remarkable that a 
number of ZEF that have very few or no neighbouring pulsars and galactic SNRs become 'selected' by the 
AGASA UHECRs (namely, ZEF No. 7, 12-17, 23, 25, 28, and 36) with the C2 triplet lying exactly inside 
ZEF No. 7.^ This observation motivated us to compare coordinates of the ZEF with positions of the nearby 
active galactic nuclei and interacting galaxies (IGs). 

4.5. Interacting Galaxies and Active Galactic Nuclei. It has been demonstrated in [33] that a number of 
particles with energies E < 10® GeV/nucleon that have had time to arrive to our galaxy from distances 
greater than 100 Mpc is completely negligible. Thus we restrict our analysis to the sources located at 
redshift z < 0.01, which corresponds to c? « 40 Mpc assuming the Hubble constant H = 75 kmsec"^ Mpc~^. 
As it will be demonstrated below, in many cases we need far less values of z to find an AGN near a ZEF. 

Figure 7 presents positions of the AGN and interacting galaxies at redshifts z < 0.01. Totally, 60 of 119 
AGN that have redshifts z < 0.005 and are located in the region 6 > —17° belong to the 4°-neighbourhoods 
of the ZEF. Twenty of them lie inside the corresponding ZEF. The same numbers for 249 AGN at redshifts 
z < 0.01 are equal to 106 and 32 objects respectively. As for 68 IGs, the corresponding numbers equal 26 
and 10, see Table 5 for a list of some of these objects. In its turn. Table 6 presents parameters of ZEF 
'selected' by the AGN and interacting galaxies. 

Let us discuss these results in more details. First of all, as is clear from Fig. 7, a number of the ZEF 
cover regions in the vicinity of the Supergalactic plane. ZEF No. 1, 7, 10, 11, 23, 28, and 37, which do not 
have neighbouring galactic sources of the types discussed above, become 'selected' by the AGN and/or IGs. 
Further 'support' is provided to ZEF No. 2, 3, 4, 13, 15, and 17. Perhaps, the most exciting observation is 



2 A recent report by Abbasi et al. [32] has revealed that one of the UHECRs with E > 3.0 x 10^^ eV 
recorded by the HiRes experiment is located at a = 169.0° and 5 = 55.9° and thus also lies inside ZEF 
No. 7. 
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Figure 7: The ZEF, AGN (A), and interacting galaxies (o) at redshift z < 0.01 [27]. Filled triangles 
and filled circles mark respectively AGN and IGs that lie within 4°-neighbourhoods of the ZEF. 
The n-like curve shows the Supergalactic plane. 



that there are numerous AGN and IGs within and near huge but 'underfilled' zones No. 14 and 30 with a hole 
bunch of galaxies located between ZEF No. 27 and No. 30. In sum, ZEF No. 14 has 23 neighbouring AGN 
and 8 IGs with respectively 19 and 7 of them at redshifts z < 0.005. Twenty-nine AGN and 2 IGs belong to 
the 4°-neighbourhood of ZEF No. 30. Fourteen of these AGN and both IGs have z < 0.005. Finally, all three 
IGs and 8 of 11 AGN that He in the 4°-neighbourhood of ZEF No. 27 are located at redshifts z < 0.005 (with 
the M58 galaxy a,t z — 0.005047). What seems to be even more important is that almost all galaxies selected 
as being close to ZEF No. 27 belong to the Virgo cluster of galaxies. These are NGC 4762, NGC 4550, 
NGC 4639, NGC 4477, M58, M88, M89, M90, M91 with the famous M87 galaxy located at A = 4.1° from 
ZEF No. 27. Galaxies M84, NGC 4380, and NGC 4438, which also belong to the Virgo cluster, lie closer to 
ZEF No. 30. Remarkably, NGC 4380 is a BL Lac object. In addition, one can notice that the Coma cluster 
and Coma supercluster of galaxies lie inside ZEF No. 14. (Their coordinates are a ^ 195°, S « 28° and 
a « 186°, S « 24° respectively [36].) In our opinion, all this makes a basis for a conjecture that zones No. 30, 
27, 14, and possibly No. 13 may be considered as parts of one huge ZEF, which to big extent originates from 
a contribution to the overall CR flux from the Coma- Virgo cluster of galaxies. 

As we have already mentioned above, AGN and IGs provide additional 'support' to comparatively large 
ZEF No. 2 and No. 3 'selected' by pulsars but 'select' a large zone No. 37 for the first time. On the other 
hand, it is hard to believe that AGN and/or IGs or even their compact groups can produce small ZEF of 
CRs in the PeV region as those No. 4, 11, 13, 15, 17, 23, and 28. It seems more likely that these ZEF should 
originate from a contribution from galactic sources. Still, we have seen that ZEF No. 11, 23, and 28 do not 
have neighbouring objects of the types considered above. The same is true for ZEF No. 7. Besides this, one 
cannot help mentioning that ZEF No. 7, 13, 15, 17, and 28 are 'selected' by both the UHECRs registered 
with the AGASA array on the one hand and AGN and IGs (except for ZEF No. 13) on the other. Thus 
we are lead to the following conjectures: (i) the corresponding UHECRs originate from the AGN selected 
in the neighbourhoods of these ZEF, (ii) a part of UHECRs are neutral particles or there are certain voids 
in intergalactic magnetic fields, (iii) a part of ZEF, especially those far from the galactic plane, are formed 
with a contribution of extragalactic CRs. 

Notice also that another AGASA UHECR lies inside ZEF No. 14 and has a number of AGN around, cf. 
Figures 6 and 7. On the other hand, three other AGASA UHECRs lie within ZEF No. 9 and 32, which have 
numerous pulsars inside but just one AGN each, both located comparatively far from the arrival directions 
of the UHECRs. In our opinion, this fact may witness in favour of an idea that a part of UHECRs may 
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Table 5. A list of some AGN and interacting galaxies (IGs) located at redshifts z < 0.01 and lying 
within 4°-neighbourhoods of the ZEF. Types of the AGN: 'Sy' is for Seyfert galaxies, possibly with 
their type (1 or 2), 'LINER' is for Low-Ionization Nuclear Emission-line Regions, and 'BLL' is for 
a BL Lac object. For ZEF with nearby AGN, the number of nearby AGN at redshifts z < 0.005 
and z < 0.01 respectively arc given in parenthesis. Names of AGN beginning with 'J' correspond 
to the Veron-Cetty & Vcron catalogue [34]. For IGs, we give nomenclature used in the SIMBAD 
database [27] with data given according to [35]. 
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originate from pulsars located in the Galaxy. This idea has already been discussed elsewhere, see, e.g., [37] 
and references therein. 

4.6. ZEF No. 16. Finally, let us briefly discuss ZEF No. 16. Recall that this ZEF has just one distant pulsar 
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Table 6. Parameters of ZEF 'selected' by AGN and interacting galaxies at redshifts z < 0.01. See 
Table 1 for notation. 
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at its left border, see Fig. 2, and no other potential sources of CRs nearby but more than 100 WR stars 
located inside its small central cell. (Also notice that there is an AGASA UHECR located to the right of 
ZEF No. 16, see Fig. 6.) Still, the zone lies in the Supcrgalactic plane. A key to a possible explanation of 
the appearance of ZEF No. 16 comes from an observation that all these WR stars belong to the M33 galaxy, 
a prominent member of the Local Group of galaxies (a « 23.46°, ^ w 30.66°; d = 730 ± 168 kpc [38]). What 
seems to be really important is that M33 hosts the most luminous steady X-ray source in the Local Group, 
see [39] and references therein, a number of radio sources, and more than 80 HII regions [36] including 
NGC 588, NGC 595, and NGC 604, a giant diffuse nebula, nearly 1500 light-years across. Thus we expect 
that ZEF No. 16 may originate from a contribution from these objects. 

Interestingly, one of successive air shower events in the 50 TeV-10 PeV energy range reported by the 
LAAS group [40] has the averaged 'coordinates' a = 21.75°, S = 30.81° and thus also lies within the central 
cell of ZEF No. 16. 

4. 7. 'Empty' ZEF. We are thus left with two ZEF without any object of the types discussed above in their 
vicinity. These are ZEF No. 19 and 22. There naturally appears a question of why are they 'empty'? 



Table 7. Parameters of 'empty' ZEF No. 19 and 22. 
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Both 'empty' ZEF arc made of regular cells. (They join into one bigger cell if we consider irregular 
cells.) Parameters of these two ZEF are listed in Table 7. One can see that there is nothing special, perhaps 
excluding the fact that the values of deviation D are close to the lowest allowed boundary. Still, this is not 
unusual for non-empty ZEF either, as is clear from Tables 1-4 and 6. It is possibly more important that 
ZEF No. 19 intersects the boundary a = 360°, and ZEF No. 22 lies close to it. At the moment, we cannot 
completely exclude a possibility that the appearance of ZEF No. 19 and 22 is an artifact of the algorithm 
resulting from a boundary effect. It is remarkable though that two of the eleven point sources in the 100 GeV 
to 100 TeV energy range with an excess greater than 4(7, reported recently by the Milagro Collaboration [41] 
are located inside ZEF No. 19. In addition, averaged coordinates of another successive air shower event 
reported by the LAAS group [40] lie within the boundaries of ZEF No. 19. 

5. EAS with TVe > 2 • 10^ 

It is undoubtedly interesting to compare arrival directions of CRs with energies just below the knee with 
those above it. Unfortunately, as we have already explained in Sec. 2, the data set under consideration 



14 



mostly covers the interval of energies just below the knee. Still, we find it interesting to perform a brief 
analysis of arrival directions of 142,223 EAS with 9 < 45° and Nc > 2 ■ 10^, which corresponds to energies 
E>2 - 10^5 eV. 




20 40 60 80 100 120 140 160 180 200 220 240 260 280 300 320 340 360 

a" 

Figure 8: Mutual positions of the ZEF obtained for the whole data set and for EAS with A'^o > 2 • 10^ 
(rectangles in thin and thick lines respectively.) 



After applying the same procedure of data analysis (with the same random seeds used at the stage of 
the data alignment) we are left with 111,370 EAS and 21 new ZEF combined of 147 regular cells, see Fig. 8. 
We do not think that 1.1 • 10^ showers are sufficient to make definite conclusions on a distribution of their 
arrival directions and will thus only give brief comments. 

First of all, one can see that there is a huge new ZEF that partially overlaps with 'old' ZEF No. 14, 27, 
and 30. This new ZEF includes all AGN and IGs selected earlier that belong to the Virgo cluster of galaxies 
and those in the region of the Coma cluster. It seems to be quite natural since one may expect that AGN 
normally produce more energetic particles than pulsars and SNRs. Next, there is another huge ZEF which 
includes 'old' ZEF No. 25 and a part of ZEF No. 26. This ZEF occurs to be mostly filled with OSCs and 
pulsars. One can also see a ZEF that is approximately adjacent to ZEF No. 6 and 9. Both SNRs located 
near the top left corner of ZEF No. 6 (G126. 2+1.6 and G127.1+0.5) belong to this new ZEF. Besides this, 
it also contains the famous Tycho SNR (SN1572). There are also a number of other new ZEF that lie close 
to the old ones. Finally, one can see a number of new ZEF that occupy areas of their own. All new ZEF 
have these or that astrophysical objects of the types discussed above in their vicinity. 



6. Conclusions 

We have presented the results of an analysis of arrival directions of more than 10^ EAS in the energy 
range « 0.1-10 PeV registered with the EAS-1000 Prototype Array. An iterative algorithm developed for 
this investigation allowed us to find a number of zones with an excess of CR flux in the data set under 
consideration at > 3a level. The zones are mostly located in the vicinity of the galactic and Supergalactic 
planes. The majority of the ZEF have powerful galactic sources of CRs inside or nearby, namely, SNRs 
(with the Crab Nebula, Cassiopeia A, and the Monogem Ring among them), pulsars, open star clusters, OB- 
associations, and regions of ionized hydrogen. This observation seems to witness in favour of the traditional 
point of view that objects of all these types make a contribution to the flux of PeV cosmic rays. 

On the other hand, there are a number of ZEF that have none or just a few objects of the above types in 
their vicinity. We find that one of such ZEF contains the M33 galaxy inside, while the majority of others have 
neighbouring active galactic nuclei and/or interacting galaxies at redshifts z < 0.01 with a big group of them 
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located in the Virgo cluster of galaxies. These observations together with a number of coincidences between 
arrival directions of UHECRs with E^A - 10^^ eV registered with the AGASA array and positions of some 
ZEF lead us to the following conjectures: (i) the flux of PcV cosmic rays is a superposition of two components, 
a galactic and an extragalactic ones (cf. [11]), and at least a part of the extragalactic component 'remembers' 
its origin; (ii) cosmic rays in the PeV (and possibly higher) energy range originate from astrophysical objects 
of various types; (iii) there are astrophysical sources that emit/accelerate particles in a wide range of energies 
so that a part of both PeV and EeV cosmic rays may have common origin. (This is not a new idea. One can 
find a number of models which attempt to explain the origin of CRs of very different energies with the help 
of a unique type of sources and/or mechanism, see, e.g., [28], [42], [9], and [10].) It may be rather difficult 
though to explain how extragalactic sources can form a noticeable flux of PeV CRs from directions close to 
their origin. Perhaps, the most obvious ideas are that at least a part of the extragalactic component consists 
of neutral particles or that there arc certain voids in intcrgalactic magnetic fields. 

The presented results pose a number of tasks for the future investigation. First of all, it is desirable to 
perform a comprehensive analysis of the algorithm with the use of a Monte-Carlo simulation. This may help 
one to figTirc out what kind of excess of the CR flux is detected and which is not. In particular, this may 
help to explain why an excess of cosmic ray flux is detected in direction to certain possible CR sources while 
to others it is not. Next, it is desirable to study arrival directions of CRs just above the knee with greater 
statistics. To implement this task, we plan to use the data obtained in late 1980's with the EAS MSU Array 
and to employ another data set already obtained with the EAS~1000 Prototype Array. Finally, we also plan 
to study another interesting feature of the CR arrival directions revealed during the investigation. These are 
zones with 'depressed' flux, i.e., zones with a number of EAS much less than the average. We think that an 
analysis of such 'sinks' can provide some information useful for understanding the origin of cosmic rays in 
the energy range near the knee. 

Acknowledgment s 

We thank A. A. Chilingarian, A. D. Erlykin, and E. Parizot for useful and stimulating communications. The 
whole investigation would have been impossible without our colleagues who obtained the experimental data. 
MZ thanks A. Golubin, A. Habib-ur-Rahman-Khan, V. Lebedev, I. Sinelobov, and R. Zhukov for the nice 
PC, which has done the main work, and B. and C. Woodworth for the enlightening discussion. 

This research has made intensive use of the NASA/IPAC Extragalactic Database (NED), which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the Na- 
tional Aeronautics and Space Administration, resources provided by SEDS [43], and especially the SIMBAD 
database, operated at CDS, Strasbourg, France. 

Only free, open source software has been used for the investigation. 



16 



References 
[1] V.F. Hess, Phys. Z. 13 (1912) 1084. 

[2] A. Chilingarian ct al., "Detection of the high-energy cosmic rays from the Monogem Ring," Astrophys. J. 
597 (2003) L129. 

[3] T. Antoni et al. for the KASCADE Collaboration, "Large scale cosmic-ray anisotropy with KASCADE," 
Astrophys. J. 604 (2004) 687. Available from arXiv: astro-ph/0312375. 

K.-H. Kampert et al., "Cosmic rays in the 'knee'-region. Recent results from KASCADE," Acta Phys. 

Polon. B35 (2004) 1799. Available from arXiv: astro-ph/0405608. 
[4] G. Benko et al., "Search for Sources of Primary Cosmic Rays at Energies Above 0.1 Pev at Tien Shan," 

Izv. RAN, Ser. Fiz. 69 (2004) 1599. Available from arXiv: astro-ph/0502065. 
[5] V.A. Kozyarivsky et al., "Mean diurnal variations of cosmic ray intensity as measured by Andyrchi Air 

Shower Array," arXiv: astro-ph/0406059. 
[6] R.W. Clay et al., "Anisotropics between 10^'^ eV and lO^* eV," Proc. 25th ICRC, Durban, 1997, 4, 185. 
[7] G.V. Kulikov, G.B. Khristiansen, Zh. Eksp. Teor. Fiz. 35 (1958) 635. 

[8] E. Roulet, "Astroparticle Theory: Some New Insights into High Energy Cosmic Rays," Int. J. Mod. 

Phys. A19 (2004) 1133. Available from arXiv: astro-ph/0310367. 
[9] S.D. Wick, CD. Dermer, A. Atoyan, "High-Energy Cosmic Rays from Gamma-Ray Bursts," Astropart. 
Phys. 21 (2004) 125. Available from arXiv: astro-ph/0310667. 
[10] A. Dar, "The Origin of Cosmic Rays— A 96- Year-Old Puzzle Solved?," arXiv: astro-ph/0408310. 
[11] H. Muraishi, S. Yanagita, T. Yoshida, "Galactic modulation of extragalactic cosmic rays: Possible origin 

of the knee in the cosmic ray spectrum," arXiv: astro-ph/0502132 ; v. 2. 
[12] T.K. Gaisser, "Origin of cosmic radiation," AIP Conf. Proc. 558 (2001) 27. Available from arXiv: 
astro-ph/0011524. 

[13] A.D. Erlykin, A.W. Wolfendale, "A single source of cosmic rays in the range 10^^-10^^ eV," J. Phys. G: 
Nucl. Part. Phys. 23 (1997) 979. 

— , "Structure in the cosmic ray spectrum: an update," J. Phys. G: Nucl. Part. Phys. 27 (2001) 1005. 
[14] S.E. Thorsett et al., "Pulsar PSR B0656-hl4, the Monogem Ring, and the origin of the 'knee' in the 

primary cosmic ray spectrum," Astrophys. J. 592 (2003) L71. Available from arXiv: astro-ph/0306462. 
[15] N. Hayashida et al., "Updated AGASA event list above 4 x 10^^ eV," Astrophys. J. 522 (1999) 225. 

Available from arXiv: astro-ph/0008102. 
[16] Yu.A. Fomin et al., "New results of the EAS-1000 Prototype operation," Proc. 26th ICRC, Salt Lake 

City, 1999, 1, 286. 

[17] L. Sun, S. Sun, "The cosmic ray incident directions observed from the Northern and Southern hemi- 
spheres," Proc. 25th ICRC, Durban, 1997, 4, 165. 

[18] J.W. Eaton, GNU Octave: A High-Level Interactive Language for Numerical Computations, Edition 3 
for version 2.0.13, 1997. (http://www.octave.org/) 

[19] D.A. Green, "A Catalogue of Galactic Supernova Remnants (2004 January version)," MuUard Radio 
Astronomy Observatory, Cavendish Laboratory, Cambridge, United Kingdom 

(available at http://www.mrao.cam.ac.uk/surveys/snrs/); see also D. A. Green, "Galactic Super- 
nova Remnants: an Updated Catalogue and Some Statistics," Bull. Astron. Soc. India 32 (2004) 335. 
Available from arXiv: astro-ph/0411083. 

[20] H.J Volk, "TeV gamma-ray observations and the origin of cosmic rays III," arXiv: astro-ph/0312585. 

[21] E.G. Berezhko, H.J Volk, "Direct evidence of efhcient cosmic ray acceleration and magnetic field am- 
plification in Cassiopeia A," arXiv: astro-ph/0404203. 

[22] T. Kobayashi et al., "High energy cosmic-ray electrons beyond 100 GeV," Proc. 26th ICRC, Salt Lake 
City 1999, 3, 61. 

K. Yoshida et al., "The origin of high energy cosmic-ray electrons and nearby supernova remnants," 
Proc. 28th ICRC, Tsukuba, 2003, 1, 1993. 
[23] A.D. Erlykin, A.W. Wolfendale, "The origin of the knee in the cosmic-ray energy spectrum," arXiv: 
astro-ph/0103477. 

[24] The ATNF Pulsar Database, http://www.atnf.csiro.au/research/pulsar/psrcat/; see also: 

R.N. Manchester, G.B. Hobbs, A. Teoh, M. Hobbs, "The ATNF Pulsar Catalogue," arXiv: astro- 
ph/0412641. 



17 



[25] J.H. Taylor, J.M. Cordes, Astrophys. J. 411 (1993) 674. 

[26] A.D. Erlykin, A.W. Wolfendalc, "Cosmic rays and the Monogem supernova remnant," Astropart. Phys. 

22 (2004) 47. Available from arXiv: astro-ph/0404530. 
[27] The SIMBAD database, http://siinbad.u-strasbg.fr/Siinbad 

[28] E. Parizot, "Superbubbles & the Galactic evolution of **Li, Be and B," Astron. Astrophys. 362 (2000) 
786. 

E. Parizot, A. Marcowith, E. van der Swaluw, A.M. Bykov, V. Tatischeff, "Superbubbles and Energetic 

Particles in the Galaxy. I: Collective effects of particle acceleration," arXiv: astro-ph/0405531. 

E. Parizot, "Cosmic-rays: an unsolved mystery at all energies!," arXiv: astro-ph/0501274. 
[29] R. Paladini et al., "A radio catalog of Galactic HII regions for applications from decimeter to millimeter 

wavelengths," Astron. & Astrophys. 397 (2003) 213. 
[30] O.V. Vedeneev et al., "EAS clusters with A^e ~ 5 x 10^," Izv. Ros. Akad. Nauk, Ser. Fiz. 65 (2001) 

1674. 

[31] Yu.A. Fomin et al., "Clusters of EAS with electron number ^ 10^," Preprint SINP MSU 2002-9/693. 

Available from arXiv: astro-ph/0203478. 
[32] R.U. Abbasi, T. Abu-Zayyad, J.F. Amann et al., "Search for Point Sources of Ultra-High Energy Cosmic 

Rays Above 40 EeV Using a Maximum Likelihood Ratio Test," arXiv: astro-ph/0412617. 
[33] G. Tarlc et al., "Limitations of the transport of cosmic ray antimatter from distant galaxies," Proc. 25th 

ICRC, Durban, 1997, 4, 205. 
[34] M.P. Veron-Cetty, P. Veron, "Quasars and Active Galactic Nuclei," Astron. & Astrophys. 412 (2003) 

399. 

[35] H. Arp, "Atlas of peculiar galaxies," Astrophys. J., Suppl. Ser. 14 (1966) 1. 

[36] The NASA/IPAC Extragalactic Database, http://nedwww.ipac.caltech.edu/index.html 

[37] M. Takeda et al., "Small-scale anisotropy of cosmic rays above 10^^ eV observed with the Akeno Giant 

Air Shower Array," Astrophys. J. 522 (1999) 225. Available from arXiv: astro-ph/9902239. 
[38] A. Brunthaler, M.J. Reid, H. Falcke, L.J. Greenhill, C. Henkel, "The Geometric Distance and Proper 

Motion of the Triangulum Galaxy (M33)," Science 307 (2005) 1440. Available from arXiv: astro- 

ph/0503058. 

[39] G. Dubus et al., "High resolution Chandra X-ray imaging of the nucleus of M33," arXiv: astro- 
ph/0406310. 

[40] N. Ochi, T. Wada, Y. Yamashita et al., "Anisotropy of Successive Air Showers," Nucl. Phys. B (Proc. 
Suppl.) 97 (2001) 173. 

[41] P.M. Saz Parkinson for the Milagro Collaboration, "Detection of TeV gamma-rays from extended sources 

with Milagro," arXiv: astro-ph/0503244. 
[42] A.D. Erlykin, T. Wibig, A.W. Wolfendalc, "A universal origin for cosmic rays above lO'^ GeV?," New 

Journal of Physics 3 (2001) 18.1. (http://www.njp.org/) 
[43] Students for the Exploration and Development of Space (SEDS), available at http://www.seds.org/. 



18 



